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ABSTRACT 

The entropy of the intracluster medium at large radii has been shown recently to 
deviate from the self-similar scaling with temperature. Using TV-body /hydrodynamic 
simulations of the ACDM cosmology, we demonstrate that this deviation is evidence 
that feedback processes are important in generating excess entropy in clusters. While 
radiative cooling increases the entropy of intracluster gas, resulting in a good match 
to the data in the centres of clusters, it produces an entropy-temperature relation 
closer to the self-similar scaling at larger radii. A model that includes feedback from 
galaxies, however, not only stabilises the cooling rate in the simulation, but is capable 
of reproducing the observed scaling behaviour both in cluster cores and at large radii. 
Feedback modifies the entropy distribution in clusters due to its increasing ability at 
expelling gas from haloes with decreasing mass. The strength of the feedback required, 
as suggested from our simulations, is consistent with supernova energetics, providing 
a large fraction of the energy reaches low-density regions and is originally contained 
within a small mass of gas. 
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1 INTRODUCTION 

Entropy has become the standard quantity for describing 
the distribution of the intracluster medium. In a self-similar 
population of clusters, where entropy originates solely due 
to gravitational infall via shock-heating, entropy scales pro- 
portional to the system's virial temperature, S oc T v i r . Such 
a model, however, does not reproduce the observed X-ray 
properties of groups and clusters. The most striking example 
of this failure is the X-ray luminosity-temperature (Lx — Tx) 
relation: the self-similar scaling, Lx oc Tx, where a = 2 
(Kaiser 1986) is flatter than the observed relation for clus- 
ters, a ~ 3 (e.g. Edge & Stewart 1991; Mushotzky & Scharf 
1997). 

This deficit in luminosity, primarily in low-mass sys- 
tems, is due to an excess of entropy in cluster cores (Evrard 
& Henry 1991; Kaiser 1991), a phenomenon first measured 
by Ponman, Cannon & Navarro (1999). Understanding the 
origin of this entropy is central to our understanding of clus- 
ter physics (Bower 1997; Voit et al. 2002, 2003). Directly 
heating the intracluster gas is the most common solution to 
this problem, however a variety of theoretical studies con- 
clude that the obvious candidates - supernovae - are at best 
only marginally capable of generating the required excess en- 



* E-mail: s.t.kay@sussex.ac.uk 

1 We define entropy as S = fcT/rig -1 keV cm 2 , where 7 = 5/3 is 
the ratio of specific heats for a monatomic ideal gas. 



tropy (e.g. Balogh, Babul & Patton 1999; Kravtsov & Yepes 
2000; Loewenstein 2000; Wu, Fabian & Nulsen 2000; Bower 
et al. 2001; Borgani et al. 2002). More energetic forms of 
supernovae (so-called hypernovae) may be important. Al- 
ternative sources of energy that are clearly observed to be 
interacting with the intracluster medium are Active Galactic 
Nuclei (AGN), capable of distributed heating through their 
release of bubbles and jets (e.g. Quilis, Bower & Balogh 2001; 
Briiggen & Kaiser 2002; Omma et al. 2003). 

Recently, attention has also focused on the effects of 
radiative cooling, which selectively removes the low-entropy 
material to form galaxies, causing higher-entropy material to 
flow in to replace it (Knight & Ponman 1997; Pearce et al. 
2000; Bryan 2000; Muanwong et al. 2001, 2002; Dave, Katz 
& Weinberg 2002; Wu & Xue 2002). Cooling is an appealing 
mechanism as it involves no free parameters (other than the 
metallicity of the gas, which is observable): the observed 
level of excess entropy is consistent with the removal of low- 
entropy gas with cooling times shorter than the age of the 
Universe, scaling as S oc T 2/3 (Voit & Bryan 2001). Cooling, 
if unchecked however, overproduces the mass in galaxies, 
requiring a feedback mechanism in order to regulate itself 
(White & Rees 1978; Cole 1991; Balogh et al. 2001), but 
as was pointed out by Voit & Bryan (2001), their argument 
still holds so long as feedback efficiently transports cold gas 
away from cluster cores. 

Perhaps the likely source of excess entropy is from a 
combination of cooling and heating processes (Kay, Thomas 
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& Theuns 2003; Tornatore et al. 2003). For example, Kay et 
al. (2003) studied the combined effects of cooling and feed- 
back (by reheating cold galactic gas) in simulations of groups 
and concluded that feedback, as well as regulating cooling, 
could contribute to the excess core entropy in systems with 
virial temperatures smaller than the heating temperature. 
In larger systems where this temperature inequality was re- 
versed, the inability of reheated gas to escape from the halo 
reduces the average entropy of the gas, and if the heating 
temperature is low enough, will affect the core entropy. 

Attention has now shifted to observing the entropy of 
intracluster gas at larger radii. In a recent paper, Pon- 
man, Sanderson & Finoguenov (2003) measured the entropy- 
temperature relation at a significant fraction of the virial 
radius (7?soo) for a sample of 66 galaxies, groups and clus- 
ters: the Birmingham-CfA Cluster Scaling Project (see also 
Sanderson et al. 2003; Sanderson & Ponman 2003). In their 
study, they measured the same deviation from self-similarity 
(S oc T 2 / 3 ) already observed in the core. This scaling be- 
haviour was also verified recently by Pratt & Arnaud (2003) , 
using XMM-Newton observations of Abell 1413 & 1983. In 
this paper, we will argue that these observations are evi- 
dence that feedback plays a significant role in forming the 
entropy distribution in clusters. 

The rest of this paper is organized as follows. In Sec- 
tion 2 we describe the simulations used to perform this 
study. Our results are presented in Section 3 and summa- 
rized in Section 4. 



2 METHOD 

2.1 Simulation details 

We have performed four simulations of the ACDM cosmol- 
ogy, setting n m = 0.3, fi A = 0.7, fl b = 0.045, h = 0.7 
& as — 0.9, consistent with the WMAP results. Each run 
started from the same initial conditions: a regular cubic 
mesh of 2,097,152 (128 3 ) particles (each of gas and dark 
matter) within a box of comoving length 60 /i -1 Mpc. This 
set the gas and dark matter particle masses to be m gas ~ 
1.3 x 10 9 and mcrai ~ 7.3 x 10 9 h' 1 M Q respectively. The 
runs were started at z = 49 and evolved to z = 0, using a 
gravitational softening length equal to 20 ft -1 kpc in comov- 
ing co-ordinates. 

The code used in this study is a parallel (MPI) ver- 
sion of GADGET (Springel, Yoshida & White 2001), an N- 
body/hydrodynamics code that uses a PM-tree to calcu- 
late gravitational forces and Smooth Particle Hydrodynam- 
ics (SPH, e.g. Monaghan 1992) to calculate gas forces. This 
version of GADGET uses entropy as the state variable in the 
time-integration of the gas (Springel & Hernquist 2002). 

Differences between each of the four runs are solely due 
to the processes incorporated that are able to change the 
entropy of the gas. These details are described below. 



2.1.1 Non-radiative model 

In this run, gas was only able to increase its entropy through 
shocks. This model reproduces self-similar scalings rather 
well. 



2.1.2 Radiative model 

Gas in this model could also lose entropy through radia- 
tive processes. We implemented the isochoric approxima- 
tion described in Thomas & Couchman (1992), using equi- 
librium cooling tables from Sutherland & Dopita (1993). We 
fixed the metallicity of the gas to Z = 0.3Z@. Dense gas 
(n > 10~ 3 cm~ 3 ,<5 > 100) which cooled below T = 1.2xl0 4 K 
was converted into collisionless 'stars'. This model signif- 
icantly overproduces the mass fraction of cooled baryons 
(observed to be 5-10 per cent; Cole et al. 2001; Balogh et al. 
2001), producing a global cooled fraction of ~ 34 per cent 
at z = 0. 

2.1.3 Feedback model 

To stabilize the cooling rate, we only allow some of the 
cooled gas to form stars and reheat the rest of the cold ma- 
terial. Similarly to Kay et al. (2003) we incorporate feedback 
effects on a probabilistic basis. (A more sophisticated heat- 
ing mechanism is unwarranted for simulations of this reso- 
lution, however this method is designed to capture the gross 
behaviour of gas being transported out of galaxies.) Rather 
than instantaneously convert all cooled particles into stars, 
we assign a value, /heat to each cooled gas particle, defined 
to be the fractional mass of cooled material that is reheated, 
and draw a random number, r, uniformly from the unit in- 
terval. If r < /heat then the particle is heated, otherwise 
it becomes collisionless. Rather than heating particles to a 
fixed temperature, we instead assigned them a fixed entropy, 
as it is this quantity that determines how far the gas can rise 
buoyantly out of the cluster. 

We ran two simulations with feedback, hereafter re- 
ferred to as Weak Feedback and Strong Feedback runs. In the 
former case we set /heat = 0.5 (i.e. a cooled particle has an 
equal chance of forming stars or being reheated) and supply 
an entropy Sheat = 100 keV cm 2 to each reheated particle. 
For the Strong Feedback simulation, we set /heat =0.1 and 
Sheat = 1000 keV cm 2 . For a cold gas particle at the same 
density, both runs require comparable amounts of energy 
per mass of collisionless material but in the Strong Feedback 
case, the energy is added to one fifth of the mass of gas in 
the Weak Feedback case. 

Note that an associated temperature depends on the 
density of the gas. At n — 10 -3 cm -3 (where most of the gas 
is heated in the simulation) both runs require ~ 1 keV of en- 
ergy per particle, comparable with the energy available from 
supernovae. However, cold gas within galaxies (our simula- 
tions do not resolve the internal structure of these objects) is 
generally at much higher densities, where much higher tem- 
peratures would be required to reach Sheat (if the gas flows 
out adiabatically) . The solution to this problem is either 
that a large fraction of energy is transported to low-density 
regions in kinetic form (e.g. Strickland & Stevens 2000) or 
that more energetic phenomena are at work, for example 
hypernovae or AGN or both. 

2.2 Cluster selection 

Clusters in each of the simulations were selected using the 
method described in detail in Muanwong et al. (2002). A 
minimal-spanning tree is first created, of all dark matter 
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Figure 1. Entropy at fixed overdensity versus temperature for 
clusters in the Non-radiative (open symbols) and Radiative (filled 
symbols) runs at Z = 0. Squares represent entropy values at 
O.IR200 an d triangles at Rsoo- Solid lines represent least-squares 
fits to the simulated data at -Rsoo- The data-points are results 
from Ponman et al. (2003), with the dashed line illustrating the 
self-similar slope (S oc T), normalized to their 8 hottest clusters. 
The horizontal dotted line is the entropy floor, S ~ f 24 keV cm 2 , 
from Lloyd-Davies et al. (2000). 



particles whose overdensity exceeds the virial value S ~ 
178fi" a55 (z) (Eke, Navarro & Frenk 1998). The tree is then 
pruned using a linking length, / = O^S" 1 ^ 3 times the mean 
interparticle separation. Spheres are then grown around the 
position of maximum density in each clump of particles un- 
til the enclosed mean density (of all particles) exceeds A 
times the comoving critical density. In this paper we focus 
on results for A = 500 and A = 200. Finally, clusters with 
fewer than 500 particles each of baryons and dark matter 
are discarded. 



3 RESULTS 

3.1 Entropy- Temperature relation 

We begin by showing results from our Non-radiative and 
Radiative simulations. Fig. illustrates entropy at a fixed 
density contrast against soft-band X-ray temperature for 
clusters in these two runs. The temperatures are calculated 



Tx = 



'E i mipiA Eoit (Ti, Z)Tj 
T.imipiA BO { t (Ti, Z) 



(1) 



where m%, pi and Ti are the mass, density and temperature of 
hot gas particle i (T< > 10 5 K) and Z = 0.3Z Q . We adopt the 
soft-band cooling function, A so f t , from Raymond & Smith 
(1977) for an energy range 0.3-1.5 keV. We calculate a mass- 
weighted entropy for each cluster by averaging the entropy 
of hot gas particles within a shell of width 20 h" 1 kpc (our 
results are insensitive to changing the size of this shell by a 
factor of 2), centred on the radius of interest. 

Squares in the figure represent entropy values at 
O.I-R200 , a radius typically chosen to highlight excess entropy 




Figure 2. Entropy at fixed overdensity versus temperature for 
clusters in the Weak Feedback (top panel) and Strong Feedback 
(bottom panel) runs at z = 0. Squares represent entropy values at 
O.I-R200 and triangles at -Rsoo- Solid lines represent least-squares 
fits to the simulated data. The data-points are results from Pon- 
man et al. (2003), with the dashed line illustrating the self-similar 
slope (S a T), normalized to their 8 hottest clusters. The hori- 
zontal dotted line is the entropy floor, S ~ 124 keV cm 2 , from 
Lloyd-Davies et al. (2000). 



in groups (Ponman et al. 1999) . As shown previously (Muan- 
wong et al. 2002; Kay et al. 2003), the Non-radiative model 
(open symbols) agrees well with the self-similar prediction 
(S oc T), apart from at low temperatures where there is 
some scatter due to insufficient resolution in the core, while 
the Radiative model (filled symbols) is in reasonable agree- 
ment with the observations (the error-bars shown here are 
from Ponman et al. 2003). 

We also show in this figure (triangles) results for a 
lower density contrast of 500, corresponding to ~ 2/3i?200- 
Again, our Non-radiative results agree well with the self- 
similar scaling but (allowing extrapolation) are around 40 
per cent higher than the expected result when normalized 
to the hottest 8 clusters in the Ponman et al. sample (up- 
per dashed line) . Note however, that the Radiative results at 
-R500 are higher than the observational data at all temper- 
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atures, but scale almost self-similarly (S oc Tx' 9 )- Cooling 
increases the entropy in clusters at all radii but in this case, 
disagrees with the observed slope. It is therefore apparent 
that cooling has been too efficient at generating excess en- 
tropy in the most massive clusters. 

Fig. |21 illustrates the same relationships for our Weak 
Feedback and Strong Feedback runs. In the former case, we 
see that this model does a reasonable job at -R500, flatter 
than the self-similar relation, but the entropy at 0.1i?2oo is 
considerably lower than observed. In this model, the gas is 
not receiving enough entropy to escape from the inner (X- 
ray dominant) region of the halo, but the feedback is having 
a larger effect on the gas in lower-mass systems, even at 
-R500. (Note that reheated gas can lose some of its entropy 
through cooling.) 

The Strong Feedback model, however, matches both re- 
lations very well. The increase in entropy allows the gas to 
escape to larger radii than previously, but the fate of this 
gas depends on the system size. In low-temperature systems, 
a significant fraction of hot gas gains enough entropy to es- 
cape the system entirely (Kay et al. 2003), resulting in an 
entropy excess that is shown here to be present even at large 
radii. This effect diminishes in higher-temperature systems, 
where the entropy level of the gas is sufficiently high that 
less material was able to escape. 

3.2 Profiles 

To illustrate more clearly the effects of cooling and heating 
on the entropy distribution of our clusters, we plot in Fig. |3] 
entropy profiles for two systems with kTx ~ 1 & 3 keV 
respectively. At large radii in both cases, the Non-Radiative 
model has an entropy profile in reasonable agreement with 
the relation, S oc R 1 ' 1 , as expected from spherical accretion 
shocks (Tozzi & Norman 2001). Both profiles flatten within 
~10 per cent of ifeoo- 

Cooling increases the entropy of the gas at all radii (by 
removing a significant fraction of gas over the age of the uni- 
verse) except for in the very centre where the cooling time 
is very short, causing the entropy to drop sharply. These 
effects are more prominent in the smaller system, reflect- 
ing the increasing efficiency of cooling with decreasing halo 
mass. 

The Weak Feedback profiles are similar to the Radiative 
profiles, except that they are lower in normalization. The 
feedback in this case is enough to prevent some gas from 
cooling but not to move it around in the cluster significantly. 
The Strong feedback model on the other hand shows greater 
differences. In the 1 keV system the entropy profile is above 
the Radiative profile out to ~ 0.4i?200, where they both be- 
come comparable. Feedback therefore has increased the core 
entropy over and above what is possible from cooling alone, 
by moving some of the hot gas (which did not cool in the Ra- 
diative model) to larger radii. In the hotter system, however, 
the Strong Feedback profile is slightly lower than the Radia- 
tive profile out to ~ 0.2i?2oo, where it decreases further, and 
at -R500 is only slightly higher than the Non-radiative profile. 
In this system, a larger fraction of gas was retained by the 
halo due to its deeper gravitational potential well. 

We note that the Strong Feedback profiles are a bit flat- 
ter than R at R500 (supported by observations). Gen- 
erally, systems above 1 keV have slopes between 0.8-0.9. 
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Figure 3. Entropy profiles for two haloes (top panel: feTx ~ 3 
keV and bottom panel: fcTx ~ 1 keV) at Z. = 0, from the Non- 
radiative simulation (solid curve), Radiative simulation (dashed 
curve), Weak Feedback simulation (dot-dashed curve) and the 
Strong Feedback simulation (triple-dot-dashed curve). The ver- 
tical dotted line marks the softening radius and the solid vertical 
line -R500- The solid diagonal line illustrates the predicted slope 
from spherical accretion shock models, S oc R 11 . 



Whether this discrepancy is a problem with this model re- 
quires further investigation, but will greatly benefit from 
both larger simulations and a larger sample of high-quality 
cluster data. 



4 SUMMARY 

In this paper we used results from iV-body/hydrodynamic 
simulations of the ACDM cosmology with various degrees 
of cooling and feedback, in order to understand the re- 
cent observational claim (Ponman et al. 2003; Pratt & Ar- 
naud 2003) that the entropy of intracluster gas does not 
scale self-similarly with temperature out to at least -R500 
(~ 2/3i?20o)- In particular, we compared our results to the 
entropy-temperature relation found by Ponman et al. at 
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O.I-R200 an d at R500, where both are found to be flatter 
(S oc T~ 2/3 ) than the self-similar scaling (S oc T). Our re- 
sults can be summarized as follows. 

• Non-radiative clusters scale self-similarly in the 
entropy-temperature plane to good approximation, both at 
O.I-R200 an d a * -^500- All but the largest clusters are observed 
to have an excess of entropy relative to the self-similar model 
at both radii. 

• Radiative cooling raises the entropy in clusters at all 
radii. While this reproduces the entropy-temperature rela- 
tion well at O.I-R200, it produces a relation that is too steep 
at 7?5oo- In this case, the relation is almost self-similar, sug- 
gesting that clusters above a few keV have too much entropy 
compared to observations. This is due to the hot gas mass 
being too low at all radii. 

• A feedback model in which gas receives a sufficient 
(1000 keV cm 2 ) amount of entropy results in excellent agree- 
ment with the observations at both radii. The reheated gas 
has enough entropy to escape the X-ray core but the fraction 
of material unable to leave the cluster (particularly at large 
radii) increases with system size, flattening the entropy- 
temperature relation. In detail, the entropy profiles are a bit 
flatter (at -R500) than suggested by observations, prompting 
future investigation into whether this discrepancy is signifi- 
cant. 

We conclude that the entropy level required to match 
the observations is high, but supernovae are energetically 
capable, providing the energy is efficiently transported into 
lower density regions and contained within a small mass of 
gas. Alternatively, more energetic phenomena (hypernovae 
or Active Galactic Nuclei) may play a part. 
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